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Abstract 

The low-temperature behavior of the static magnetic susceptibihty x(^) of exchange-disordered 
antiferromagnetic spin chains is investigated. It is shown that for a relatively small and even number 
of spins in the chain, two exchange distributions which are expected to occur in nanochains of P 
donors in silicon lead to qualitatively distinct behaviors of the low-temperature susceptibility. As a 
consequence, magnetic measurements might be useful to characterize whether a given sample meets 
the requirements compatible with Kane's original proposal for the exchange gates in a silicon-based 
quantum computer hardware. We also explore the dependence of xi^) on the number of spins in 
the chain as it increases towards the thermodynamic limit, where any degree or distribution of 
disorder leads to the same low-temperature scaling behavior. We identify a crossover regime where 
the two distributions of disorder may not be clearly differentiated, but the characteristic scaling of 
the thermodynamic limit has not yet been reached. 

PACS numbers: 75.10.Jm, 71.55.Cii 
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I. INTRODUCTION 



Low dimensional quantum antiferromagnets have been intensively studied over the years, 
and became a reference problem in Condensed Matter Physics, Strongly Correlated Sys- 
tems, and Statistical Mechanics.— Quantum spin chains may exhibit quantum phase transi- 
tions and, in the continuum limit, turned out to be an active area for application of quan- 
tum field theory in Condensed Matter Physics. Random quantum antiferromagnetic (AF) 
chains, which model the magnetic behavior of several quasi-one-dimensional compounds, 
have also attracted a great deal of attentioii^. Theoretical studies^ii of the quantum spin- 
1/2 disordered Heisenberg antiferromagnet in ID show that, for infinite chains, the magnetic 
susceptibility behavior is essentially insensitive to the specific disorder distribution. While 
the behavior in the thermodynamic limit is thus well established, there has not yet been 
a detailed study of how that limit is approached, and specifically how the behavior of the 
magnetic susceptibilities of finite chains might depend on the disorder distribution. This is 
an interesting question in light of current efforts towards nanoscale control and applications 
which, coupled to improved techniques of material growth, directed interest into magnetic 
nanostructures, in particular into finite magnetic chains.— 

One of the most exciting potential applications of magnetism at the extreme micro- 
scopic level consists in utilizing the 2-level dynamics of the electron spin as a physical 
implementation of the quantum-bit (qubit) in a solid state quantum computer, where the 
required entanglement between qubits could be provided by the exchange coupling between 
electrons:^ For example, the quantum behavior of exchange-coupled electrons bound to an 
array of phosphorous donors in silicon is a key element in Kane's proposal for a Si-based 
quantum computer.'' Low-temperature magnetic susceptibility measurements of P doped 
Si have provided valuable information regarding the exchange distributions in randomly 
3D-doped samples.^ It is expected that similar measurements might also be relevant in char- 
acterizing linear arrays of donors. In Kane's proposal, two-qubit operations are mediated 
by the exchange interaction between electrons bound to nearest-neighbor P atoms in the 
chain. Previous studies have shown that this exchange coupling J is always AF, and that 
its strength is highly sensitive to the inter-donor relative positioning . Indeed, changes 
of just one lattice parameter in the relative positioning of two P impurities may alter the 
magnitude of the coupling between them by orders of magnitude. Controllable exchange 
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coupling compatible with Kane's original proposal would be achieved if all donors in the 
chain could be positioned exactly along a single [100] crystal axis.— In this situation, the 
exchange coupling behavior is (as assumed by Kane'') similar to the hydrogenic case,— i.e., 
it decays exponentially with increasing interdonor distance. Uncertainties in the interdonor 
distances along the chain would result in a narrow distribution of values for J around a "tar- 
get" value Jq. Assuming a substitutional donor positioning precision of about 1 nm in the Si 
lattice,— the perfect [100] alignment situation may be modeled by a trimodal exchange dis- 
tribution Ptri{J). However, if instead of perfect alignment the donor positions are randomly 
distributed among all substitutional sites within a small spherical region of 1 nm around the 
ideal ( "target" ) impurity sites, the peculiar band structure of Si leads to a wide distribution 
of exchange coupling, peaked at J = 0^^ causing difficulties in the operation and control of 
the "exchange gates". We approximate such distribution here by an exponential function, 
Pexp{J) for J > 0. It is therefore clear that the exchange disorder distribution is also highly 
sensitive to the positioning distribution of the P donors in Si. 

The magnetic susceptibility behavior of infinite chains is essentially insensitive to the 
specific disorder distribution, as mentioned above, indicating that susceptibility measure- 
ments would not differentiate among the two cases mentioned above in this limit. Of course 
this may be different for finite chains, a situation which is also of practical interest in terms 
of guiding current fabrication efforts towards P donors positioning in Si .-^^i-^^ The aim of 
the present work is to shed light on this problem by investigating the relation between the 
magnetic response of linear chains of spins and the distributions of the exchange interaction 
within the chains. Our results indicate that for AF disordered chains with even and rela- 
tively small number N of sites, the exchange distributions Ptri{J) and Pexp{J) lead to quite 
distinct low temperature behavior of the zero-frequency uniform magnetic susceptibility. 
Hence, the two distributions could be experimentally distinguished by sufficiently sensitive 
magnetic measurements, providing useful information regarding donor alignment. 

This paper is organized as follows. In Sec. |H] we briefly review results available in the 
literature regarding infinite AF chains. In Sec. IIIII we consider finite chains with relatively 
small number of spins, starting from spin pairs and trios for which analytical solutions are 
obtained, as well as 8-spin chains, which are solved numerically. In Sec. IIVI we analyze 
the crossover into the thermodynamic limit by solving for longer chains via a Quantum 
Monte-Carlo method. Our summary and conclusions are presented in Sec. |3 
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II. INFINITE ANTIFERROMAGNETIC CHAINS 



The Hamiltonian describing an open chain with spins is 

N-l 
i 

where the spin quantum number is S = 1/2. Since we are interested in AF chains, we 
assume that Jj > for every i. We briefly review in this section several pertinent results 
available in the literature for AF chains in the N oo limit. 

The ground state of an ordered inflnite chain {Ji = J for every i) can be obtained from 
Bethe ansatz. Grifiithsi^ has shown that the zero temperature susceptibility per spin of such 
system is flnite and given by x(T = 0)/xo{'J) = V^^i where Xoi^) = Q^^I^bYIJ, where jis 
is the Bohr magneton and g{= 2) is the Lande factor. For general T, fleld theory methods 
[k = 1 Wess-Zumino-Witten (WZW) non-linear a modelp^ give 

X„(J) (Jr)^ I 2in(r„/r) 

where Tq is a temperature cutoff. Quantum Monte-Carlo (QMC) calculations of x(T) have 
been carried out by Kim et al^ and their results for Tq/ J = 1.8 (we use units of energy for 
temperature, that is, the Boltzmann constant /c^ is set equal to one) are well fltted by the 
WZW expression ((2)). 

According to real-space renormalization group theory,— the introduction of any amount 
of disorder drives the system into a random singlet phase, in which each spin forms a singlet 
pair with another spin; pairs with arbitrarily long distance also exist. Bonds among distant 
spins, however, correspond to very weak coupling. The low-temperature excitations basically 
involve breaking these weakest bonds, resulting in nearly-free spins giving rise to a Curie 
susceptibility modifled by the statistics in the number of contributing spins. As a result the 
magnetic susceptibility at low T diverges as^ 

(3) 

III. FINITE CHAINS 

The above results refer to chains in the thermodynamic limit (A^ — ^ oo). Foi finite ordered 
rings (periodic boundary conditions), Bonner and Fisher— have calculated the susceptibility 




X(r^0)~l/ T(logT) 
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per spin Xn{T) for up to N=ll spins based on direct diagonalization of TC in the presence 
of a magnetic field. They have found that Xn{T —>■ 0) exhibits distinct behavior depending 
on whether is even or odd. In the first case, pairs of neighboring spins tend to form 
singlets and XN=even(T ^ 0) ^ 0, whereas in the second case, the occurrence of unpaired 
spins leads to a Curie-law behavior XN=oddiT ^ 0) ~ 1/T oo. These results immediately 
raise the question as to what extent the behavior of xn{T) changes by the introduction of 
disorder. Moreover, could we infer, based on the magnetic response, the type of exchange 
disorder distribution? The relevant distributions here are: (i) trimodal, with 

Pt„{J) = (l/3){5( J - Jo) + S{J - (1 + W)Jo) + S{J - (1 - W)Jo)} , (4) 

where Jo and W are both positive, with W < 1, and (ii) exponential, with 

Pe.p{J) = ^e~'/'° e(J) , (5) 

-'0 

where G is the step-function. Note that in both cases (J) = Jo, the exchange "target" value. 
We consider initially the N = 2 case, for which the susceptibility per spin is given hy^ 

X2{T, J) ^ PJ 
Xo(J) 3 + e/3^' 

where /3 = 1/T and Xo(J) is given above Eq. (j21). Considering the average of X2 over the 

oo 

above distributions, it is clear that {X2)tri = / Ptri{J) X2(T,J) dJ vanishes as T 0, as 



in the absence of disorder. For the exponential distribution, it is convenient to split the 
integral for {X2)exp ^'^^'^ terms 

{X2)exp = -r / 0J ^'^+T , Q , flj in 

Jo Jo 6 + e'^'' Jo Ja/p 6 + ef^-^ 

where a is a constant, and {x2)exp is given in units of {g^sY- For sufficiently low temper- 
atures, a can be chosen such that a/ (3 <^ Jo and << 1. Hence, in the first integral 
(corresponding to small values of J), we can approximate e~-^l^^^ 1, while in the second 
one, the term e'^"^ is always much greater than 1, leading to 

Jo(X2)e.p~Ci + e-"^Ci, (8) 

where Cx = /o° 1/(3 + e^) dx ^ C 1/(3 + e^) dx = ln(4)/3. Thus, as T 0, {x2)exp 
approaches a non-zero value, in contrast with the trimodal distribution result. 
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For an open = 3 chain, the susceptibihty per spin isi^ 

5 + e/3(Ji+J2)/2 cogh ({(3/2)Jjl-JiJ2 + ji] 

X3(T, Ji, J2) = (/5/12) ) (9) 

1 + e^(-^i+-^2)/2 cosh K/5/2)^j2- JiJ2 + ^|) 

It is clear that for non- negative values of Ji and J21 Xs exhibits a Curie-like divergence as 
T ^ 0. As a consequence, the average of xs over the trimodal distribution, {x?)trv ^^^^ 
diverges as T approaches 0. Regarding the exponential distribution, the change of variables 
Ji = J sin 6 and J2 = J cos 6 leads to 

J^AX.U. - WU) f^iof e-<.— )^/-i±W J ,J, (10) 

where f{x,e) = e^(«i'^^+=°'^^)/2 cosh ((a;/2)Vl - sin ^ cos 0). As for the A^ = 2 case, we split 
the integral above into two, which we label Ji and I2, corresponding to < J < a/ (3 and 
a/ P < J < 00, respectively. For sufficiently low T, a is again chosen such that a/13 <^ Jq 
and e"° << 1. Since 1 < sin6' + cos^ < v^, the term e^i^''-^'^+^°^'^)JIJo the first integral 
(corresponding to small values of J) can be approximated by 1, giving 

where C[ is T-independent. In the second integral, since I3J > a and a is large, the T- 
dependent term in both the numerator and denominator of the integrand is much larger 
than 1, so that 

de e-{sine+cose)j/Jo J ~ / / g-(sme+cose)j/Jo j ^ (^2) 

Jo Ja//3 Jo Jo 

where C'^ is T-independent. Thus, from Eq. Mxz)exp = {(3/12){Ii + l2)/Jo ~ I^C'J Jq, 
which diverges following a Curie law as T — 0. 

For a better physical insight, it is instructive to compare the N = 2 and A^ = 3 cases. We 
note that in the low-T regime, the T-dependent behavior of xn is dominated in the first case 
by the disorder distribution , while in the second by the odd parity of A^. The Curie-type 
low-T behavior of {xsiT)) is related to the occurrence of "unpaired" spins, independently 
of P{J)- It is interesting that the large- J tail of the distribution gives the dominant contri- 
bution to the average behavior, namely the integral given by I2 in Eq. (|12p. corresponding 
to the more strongly coupled 3-spin chains. On the other hand, for N = 2, the behavior of 
(X2(T)) at low T is dominated by the contribution from small values of J (more precisely, 
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FIG. 1: Low-temperature behavior of the magnetic susceptibihty per spin of disordered antiferro- 
magnetic spin chains with exponential and trimodal distributions of exchange coupling. The plots 
give Xn/xo{Jo) versus T/Jq for N = 2 and 8. Results for the exponential disorder distribution and 
for ordered chains are given on the left panel, whereas those for the trimodal disorder distributions 
with width W = 0.3 and 0.5 are given on the right panel. All results are for open boundary 
conditions. For N = 8 the statistical error bars are smaller than the data points, and the straight 
lines joining the points are just guides to the eye. 

from values of J <^ mm{T, Jq}), as shown in Eqs. (|7j) and The occurrence of arbitrarily 
small values of J in the exponential distribution weakens the tendency of neighboring spins 
to form singlets,— leading to a finite value of {x2(T = 0))exp, while {x2iT = 0))tri = 0. 

We expect the above considerations to apply to other finite values of A^. As odd- 
parity N chains always lead to unpaired spins, and thus to a Curie-like low-T divergence 
of {xN=odd(T)), susceptibility measurements would not be useful for distinguishing among 
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different exchange disorder distributions if is odd. Since analytical calculations become 
impractical as increases, in order to illustrate the even-parity case we have carried out 
numerical calculations for both {x8)tri and {x&)exp as a function of T. We obtain the sus- 
ceptibility numerically by determining the spin-spin correlation function (S^Sj) from which, 
through the fluctuation-dissipation theorem, we obtain x = PJ2i,j{SiSj). For each tem- 
perature we have averaged over 10,000 realizations of disordered configurations. For the 
trimodal distribution, two values for the width parameter have been considered, namely 
W = 0.3 and 0.5. Results are presented in Fig. [H which shows susceptibility curves for 
chains with N = 2 and 8, for trimodal and exponential disorder distributions, as well as for 
the ordered chain cases. We note that results for {xN{T))tri are qualitatively very similar, 
regardless of the width parameter W, even in the limit W = 0, corresponding to the ordered 
chains. All curves reach a maximum and eventually decrease as T ^ 0, going to zero for 
T = 0. We remark the obvious fact, also illustrated in Fig. ^ that as W increases the 
maximum in x and the sharp downturn toward zero value occur at lower temperatures. The 
results for the exponential distribution are markedly different, with an increasing {xN{T))exp 
for decreasing T. In principle the trimodal and exponential distributions might be identified 
and differentiated through low-temperature susceptibility measurements in such even-parity 
chains. 

IV. CROSSOVER REGIME AND THERMODYNAMIC LIMIT 

Further increase in A^ eventually leads to the thermodynamic limit behavior— given in 
Eq. (121). For even A^, approach to such behavior, giving a divergent x as T — 0, might be 
expected from the comparison between the results of {x2iT))exp and {x?,iT))exp in Fig. [H 
but it is not so clear for the trimodal distributions. Another puzzling point regards the 
sensitivity to even {Neven) or odd (Nodd) values of A^. For smaAl-Nodd, X exhibits a Curie-like 
~ 1/T divergence at T ^ independent of the disorder distribution, as discussed for A^ = 3 
in Sec. IIIH while small- A^g^en chains are quite sensitive to -P(J). But odd- and even-chains 
results must approach each other, possibly with some sensitivity to the type of disorder 
(remanent from Neven)-, as A^ increases. 

In order to clarify these points we have calculated the susceptibility for larger values of 
A^ using a stochastic series expansion (SSE) method. The SSE is a QMC method based 
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on the Taylor expansion of the Boltzmann weight operator up to a very high order.— 
The partition function and observables can then be evaluated via importance sampling of 
the different terms appearing in this series. Choosing a large enough order for the expan- 
sion, the systematic errors introduced by the truncation of the series are negligible. The 
SSE method also allows us to use importance sampling update schemes based on global 
changes of the system (cluster or loop updates) that are extremely efficient, especially for 
the highly symmetric Heisenberg model studied here. We have adjusted the precision of the 
data obtained through importance sampling so that the errors obtained on each individual 
realization are roughly one order if magnitude smaller than a typical difference between two 
realizations. We have considered open chains and averaged our results over 10,000 disorder 
realizations (except for the large size N = 128 where we used periodic boundary conditions 
and averaged over only 100 realizations). 

We present in Fig. |21 results for {xN(T))tri and {xN(T))exp for increasing up to = 17 
for the exponential and trimodal {W = 0.5) distributions. We identify here the following 
trends towards the thermodynamic limit: (i) from the frames for Ngt^en (on the left) we see 
that the well differentiated disorder distribution results for = 8 in (a) approach each other 
as Neyen iucreases [(c) and (e)]; (ii) from the frames for Nodd (on the right) we do not identify 
significative changes as the Curie behavior discussed in Sec. Illll for N=3 acquires logarithmic 
corrections [eventually leading to the asymptotic behavior governed by Eq. Q] which are not 
easily captured on this scale and through this range of Nodd values; (iii) following successive 
rows we see that the even-odd differences (namely going from Neven to Neven + 1) become 
less prominent as Neven increases [(a)-(b), (c)-(d), (e)-(f)]. From the practical point of view, 
it is clear from (i) that, contrary to the N = 8 results, for = 16 the two distributions 
lead to very similar qualitative types of behavior of the susceptibility, indicating that, for 
the particular distributions and temperature range considered here, they may not be clearly 
differentiated beyond this size of chains, regardless of the parity of A^. 

Approach to the thermodynamic limit, given in Eq. (jH)), is usually investigated by plotting 
(Tx)~^/^ vs logT/ Jo, which leads to a linear behavior at low temperatures in this limit.-^ In 
Fig. El we illustrate the crossover regime by presenting such scaled plots for A^ = 16 and 128 
for trimodal {W = 0.5) and exponential distributions. In all cases, straight lines (dashed for 
trimodal and solid for exponential) are drawn through the two lowest-T calculated points. 
We note that for the trimodal distribution and A^ = 16 such line does not include any other 
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FIG. 2: Magnetic susceptibility per spin for disordered finite chains of N spins for trimodal {W = 
0.5) and exponential disorder distributions. For T ^ the = 8 results for the two distributions 
are clearly distinct, while for = 9 they are practically collapsed. As the number of spins in the 
chain increases, the results for the two distributions and N^ven (left column) approach each other, 
as results for N^ven and N^ven + 1 chains (successive rows) also approach each other. Lines are 
guides to the eye and statistical error bars are smaller than the data points. 

calculated point, whereas for = 128 it gives a good fitting for T / Jq up to about 0.25, 
indicating that the thermodynamic limit has been reached up to this temperature for this 
value of A^. For the exponential distribution, the approach to the thermodynamic limit 
scaling with increasing A^ is faster, although careful analysis of the data plotted in Fig. El 
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FIG. 3: Scaled plots for the low-temperature behavior of the magnetic susceptibility per spin 
of disordered chains, illustrating the approach to the thermodynamic limit as N increases. The 
plots give \/Txn /xq{Jo) versus T / Jq (logarithm scale) for exponential and trimodal {W = 0.5) 
distributions in chains with = 16 and 128. The thermodynamic limit (A^ oo) corresponds 
to the linear behavior indicated by the straight solid and dashed lines drawn to connect the two 
lowest-T data points for exponential and trimodal distributions respectively. 
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shows that the points for = 16 and exponential distribution [triangles in (a)] down-turn 
from the solid line at small T, and the good agreement for T/Jq = 1 in (a) is fortuitous, 
as it does not remain for = 128 [triangles in (b)]. It is interesting to note that the solid 
(dashed) line in Fig. |21|^a) is nearly parallel to the solid (dashed) line in (b), indicating that 
some aspects of the — oo behavior are already captured at the smaller values at low-T. 
Fig. ini^b) shows that in the thermodynamic limit the susceptibility is quantitatively quite 
sensitive to the disorder distribution, although in practice this is probably not as valuable a 
tool to identify the disorder distribution as the differences encountered for the smaller even 
values of A^ [e.g. Fig. El^a)]. 

V. SUMMARY AND CONCLUSIONS 

We have investigated the low-T behavior of the uniform magnetic susceptibility of ex- 
change disordered spin-1/2 AF chains as a function of the number A^ of spins in the chain 
and disorder, for trimodal and exponential disorder distributions. Formally, the key dis- 
tinction among the two distributions considered here is that the exponential distribution is 
not bound, so that pair exchange coupling J arbitrarily close to zero may occur, while the 
trimodal distribution is bound and J does not become arbitrarily small. 

For chains with even and relatively small number of spins {Neven ~ 8) the susceptibility 
displays distinct behaviors for the two exchange distributions which are expected to occur 
in nanochains of P donors in Si. According to our results in Fig. ^ it might be possible 
to identify the atomic-scale positioning of the P atoms in such chains, thus providing com- 
plementary information on whether sample preparation techniques meet the requirements 
compatible with Kane's original proposal for a quantum computer hardware. For larger 
values of A^ (A^ ^ 16) such differentiation would not be so straightforward, as illustrated in 
Fig. El We also note that xn{T) becomes less sensitive to A^ being even or odd when A^ ^ 16. 
For completeness, we have also investigated the approach to the thermodynamic limit, and 
we have found that for A^ = 128 the expected scaling of x{T — * 0) is already obtained, while 
for A^ = 16 the system is still far from the thermodynamic limit, as discussed in Fig. El 

We therefore identify three regimes as A^ increases: (i) when A^ ^ 8 the even-A^ disor- 
dered chains present quite distinct behaviors according to the exchange distribution; (ii) for 
intermediate A^ values, illustrated here by A^ = 16, the low temperature behavior of the 
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magnetic susceptibility corresponding to different distributions is not easily differentiated, 
although the thermodynamic limit has not yet been reached; (iii) for larger A^, illustrated 
here by = 128, the two distributions follow the thermodynamic limit scaling, which is a 
signature that the random singlet phase is formed at the lowest T. 
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